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ABSTRACT 

Context. The recently developed test-field method permits to compute dynamo coefficients from global, direct numerical simulations. 
The subsequent use of these parameters in mean-field models enables us to compare self-consistent dynamo models with their mean- 
field counterparts. So far, this has been done for a simulation of rotating magnetoconvection and a simple benchmark dynamo, which 
are both (quasi-)stationary. 

Aims. It is shown that chaotically time-dependent dynamos in a low Rossby number regime may be appropriately described by 
corresponding mean-field results. Also, it is pointed out under which conditions mean-field models do not match direct numerical 
simulations. 

Methods. We solve the equations of magnetohydrodynamics (MHD) in a rotating, spherical shell in the Boussinesq approximation. 
Based on this, we compute mean-field coefficients for several models with the help of the previously developed test-field method. The 
parameterization of the mean electromotive force by these coefficients is tested against direct numerical simulations. In addition, we 
use the determined dynamo coefficients in mean-field models and compare the outcome with azimuthally averaged fields from direct 
numerical simulations. 

Results. The azimuthally and time averaged electromotive force in fast rotating dynamos is sufficiently well parameterized by the 
set of determined mean-field coefficients. In comparison to the previously considered (quasi-)stationary dynamo, the chaotic time- 
dependence leads to an improved scale separation and thus to a better agreement between direct numerical simulations and mean-field 
results. 



1. Introduction 

Mean-field electrodynamics provides conceptual understanding 
of dynamo processes generating coher ent, large-scale magnetic 
fields in planets, stars and galaxies (iKrause & Radleii 1980t 
MofFatt 1978). But, it suffers from a number of free parame- 
ters, which are only poorly known under astrophysically rele- 
vant conditions. We briefly summarize the essentials of mean- 
field theory in order to specify this point. Within a mean-field 
approach, the velocity and the magnetic field are usually split in 
large-scale mean fields, V, B, and residual paits, v, b, varying on 
much shorter length or time scales, 

V=V + v, B^B + b. (1) 
Then, the evolution of the mean field is governed by 



(|2]l only, which requires to express £ as a linear functional of 
V, V and B. This is justified with regard to Q, and we may write 



+ J J Kij(x,x',t, t')Bj(x',t')d\'dt', 



(4) 



in which K denotes some integral kernel and we refer to carte- 
sian coordinates for the moment. The term in (|4| accounts 
for small-scale dynamo action and is not considered here; it will 
be discussed later In addition, we assume that £ depends only 
instantaneously and nearly locally on B. This is a crucial as- 
sumption, it will be referred to as the assumption of scale sepa- 
ration in the following. Therefore, B in (HJl may be replaced by 
its rapidly converging Taylor series expansion at x, 



— = V X (£ -H y X B - 77V X B). 
of 



- , - , dBj(x,t) 
Bj(x', t) = Bj{x, t) + (4 - Xk) — + ■ 



(2) 



and taken out of the integral: 



In the above equation, the mean electromotive force is defined 
as £ = V xb and rj is the magnetic diffusivity. Equation ^ is 
coupled with an equation for the residual field 



- 9Bj 
£; = aijBj + bijk- — + 



y X -t- V X B - 77V X ft), 



(3) 



in which G is defined asG = vxft-vxft. The above system 
of equations is equivalent to the usual induction equation and no 
approximation has been applied so far An important simplifica- 
tion of the mean-field approach consists in considering equation 



with 



^ J Kijix, x', t, t') (f'x'df , 

J" J" Kijix, x', t, t')(x[ - Xk) d^x'dt' 



(5) 

(6) 

(7) 
(8) 
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The dynamo coefficients a and b defined in ^ and ([8]l together 
with the expansion (|6]l may finally be used to solve equation (|2]l. 
Moreover, a and b are directly linked to dynamo processes, i.e. 
to the generation, advection and diff'usion of the mean magnetic 
field, and thus provide physical concepts to explain dynamo ac- 
tion (e.g. [Radler 1995). Unfortunately, they are not known in 
general and previous work often reUes on arbitrary assumptions 

on them. 

The test-field method by 'Schrin ner etlH (l2005l |2007|) per- 
mits to compute a and b from direct numerical simulations. It 
was first applied to a simulation of rotating magnetoconvec- 
tion (Olson etal. 1999) and a simple geodynamo simulation 
(jChristensen et al. 2Q01), which are both stationary except for an 
azimuthal drift. In the past three years, the test-field method was 
intensively used to compute mean-field coefficients for box sim- 
ulations (e.g. Sur et al. 2008; Gressel et al. 2008; Brandenburg 
120091; iKapvla et alj|2009t IRadler & Brandenbu rg 2009), which 
are not in ffie focus of this paper Although ffie test-field method 
proved to be reliable, the parameterization of the electromo- 
tive force for the steady geodynamo model was not satisfactory. 
The expansion in (|6]l does not converge for th is steady exam- 
ple b ecause of a non-sufficient scale separation (ISchrinner et al.l 
l2007i) . In this paper, we revisit the problem for chaotically time- 
dependent dynamos. We test azimuthally and time averaged 
electromotive-force vectors against their parameterization based 
on corresponding dynamo coefficients and compare azimuthally 
and time-averaged fields from direct numerical simulations with 
mean-field results. 

2. Dynamo simulations 

We solve the equations of magnetohydrodynamics (MHD) in the 
Boussinesq approximation for a conducting fluid in a rotating 
spherical shell as given by Olson et al.. (.19991 

E -H y ■ VF - V^yj + 2zxV + VP = 

Ra — T + —(VxB)xB (9) 

To Pm 

^ = Vx(vxB) + — V^B (10) 
dt Pm 

— -H V ■ vr = — v^r. (11) 

dt Pr 

The coupled differential equations for the the velocity V, the 
magnetic field B, and the temperature T are governed by four 
parameters. These are the Ekman number E - v/Q.L^, the (mod- 
ified) Rayleigh number Ra = argoATL/vQ., the Prandtl num- 
ber Pr = v/k and the magnetic Prandtl number Pm - v/rj. In 
these expressions, v denotes the kinematic viscosity, Q. the ro- 
tation rate, L ffie shell width, aj the thermal expansion coeffi- 
cient, go is the gravitational acceleration at the outer boundary, 
Ar stands for the temperature difference between the spherical 
boundaries, k is the thermal and 77 = 1 /ficr the magnetic diffusiv- 
ity with the magnetic permeability fi and the electrical conduc- 
tivity cr. Besides these input parameters, we introduce the mag- 
netic Reynolds number Rm = V^ms L/rj, and the local Rossby 
number /?o/ - Vrms/ i^L) ■ (Jin) as important output parameters. 
In the latter expressions, Vnns denotes ffie rms-velocity and nfl 
is the mean half-wavelength of the flow derived from the kinetic 
energy spectrum (Christensen & Aubert 2006). 

No-slip mechanical boundary conditions were chosen for all 
simulations presented here and the magnetic field continues as a 



potential field outside the fluid shell. Convection is driven by an 
imposed temperature difference between the inner and the outer 
shell shell-boundary. 

3. Computation of dynamo coefficients 

In order to determine the dynamo coefficients a and b, we apply 
ffie tes t-field method explained in full detail in ISchrinner et al.l 
( l2007h . The principal idea of ffiis approach is to measure the 
mean electromotive force generated by the interaction of the 
flow with an arbitrarily imposed test field, Bf. The imposed test 
field appears as an inhomogeneity in equation (O, 

— -Vx(G + yxi.-77Vxfe) = VxvxBT, (12) 
dt 

which is solved simultaneously with dglfTTTi. The electromotive 
force due to a given Bj may then be computed as £t = v xb 
with b resulting from (fT2t . Finally, Bj is used to solve relation 
(|6]l for the dynamo coefficients. In order to close the resulting 
system of linear equations and to determine all components of a 
and b, the numerical experiment (fT2t must be repeated several 
times with different fields Bj. 

So far, averaged quantities labelled by an overbar are ax- 
isymmetric fields whereas v and b in ( fT2l i are non-axisymmetric 
residuals. For a stochastically stationary but nevertheless chaoti- 
cally time-dependent flow, £ and also the mean-field coefficients 
vary in time. Thus, we introduce an additional time averaging in- 
dicated by brackets, <■■■>, and write 

&^a<B> +b <VB> + ■■ (13) 

instead of (|6]l. In the above equation, £ has been time averaged, 

&^<Vxb>, (14) 

and also a and b denote time averaged tensors. Furthermore, we 
interpret averages in (|2} as combined azimuthal and time aver- 
ages. Both is formally not correct. However, the approximations 
applied can be justified, if 

i) temporal fluctuations of the axisymmetric velocity and 
magnetic field are comparatively small, and 

ii) time averages of the non-axisymmetric residuals of the 
velocity and the magnetic field are negligible. 

The above assumptions imply that azimuthal and time aver- 
aging act in a similar way and reinforce each other As a conse- 
quence, mean-field coefficients are derived in this context from 
time averages of vectors £t resulting from (fT2t . We emphasise 
that assumptions i) and ii) are not justified a priori but have to be 
tested in numerical simulations. 

4. Comparison of mean-field results with direct 
numerical simulations 

Mean-field coefficients determined from direct numerical sim- 
ulations may be used in mean-field models. Subsequently, az- 
imuthally and time-averaged magnetic fields resulting from 
three-dimensional, self-consistent models may be compared 
with their mean-field counterparts. Mean-field calculations pre- 
sented here are based on equation (|2} written as an eigenvalue 
problem, 

<tB = VxDB, (15) 
in which the linear operator D is defined as 

DB = FxB-i-aB-HbVB- 1/PmVxB, (16) 
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Table 1. Overview of the models considered, ordered with re- 
spect to their local Rossby number 



Model 


E 


Ra 


Pm 


Mean 1 


Roi 


Rm 


Afi 




1 


1 X 10"' 


100 


5 


5 


0.013 


39 


4.6 


4.18 


2 


1 X 10-* 


334 


3 


11 


0.014 


123 


1.6 


-7.05 


3 


1 X 10-* 


334 


2 


11 


0.015 


86 


1.4 


-3.87 


4 


3 X 10-* 


195 


3 


9 


0.019 


67 


1.6 


2.22 


5 


3 X 10-* 


243 


2 


9 


0.024 


56 


1.2 


0.66 


6 


3 X lO-'* 


285 


2 


9 


0.026 


61 


1.3 


0.05 


7 


3 X 10-* 


375 


1.5 


11 


0.042 


60 


1.1 


0.78 



and a and b are time-averaged dynamo coefficients as described 
above. The eigenvalues cr lead to a time evolution of each mode 
proportional to exp (erf). For more details concerning the eigen- 
value calculation, we refer to lSchrinner et al.l(l2010bl ). 

In addition, the parameterization of the electromotive force 
by dynamo coefficients, i.e. relation (|6|l, may be tested against 
the mean electromotive force in direct numerical simulations. 
Again by virtue of (11311141 1. we compare £ as given in ( fT4b and 
derived directly from numerical models with 
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Fig. 1. Volume-averaged kinetic energy density versus magnetic 
diffusion time for model 3. 



G^'> = a <B> +b <VB>, 



(17) 



in which <B> and < VB > are also obtai ned from direct numer- 
ical simulations. Furthermore, we follow ISchrinner et al.l (l2007h 
and introduce 

L \G - dV 



in order to quantify errors in the parameterization of £. 



to 
< 




5. Results 

Besides the (quasi-)steady benchmark dynamo (model 1), we 
considered 6 chaotically time-dependent dynamos. The models 
are defined by 4 control parameters; E, Ra, and Pm were varied 
(see Table[T]), whereas the Prandtl number was kept fix and set to 
1 for all simulations. The local Rossby number is always lower 
than 0.12 for the models considered here. They t herefore belong 
to the regime of kinematically stable dynamos (ISchrinner et alj 
1201 Oah . In order to simplify the time averaging, models with 
low and fairly moderate Rm were chosen. The resulting velocity 
fields are almost symmetric with respect to the equatorial plane, 
except for model 7, and convection occurs mainly outside the in- 
ner core tangent cylinder Figure [1] shows the volume-averaged 
kinetic energy density for model 3 and illustrates the irregular 
time-dependence of these models. 

The parameterization of the mean electromotive force given 
by ( [TtT i improved a lot for the chaotically time-dependent dy- 
namos in comparison to the benchmark dynamo; A£ decreases 
by more than a factor of 4 from model 1 to model 7. The consid- 
erable drop of A£ from the steady to the time-dependent models 
is most clearly visible in Fig.|2] It shows Afi for different models 
versus their magnetic Reynolds number Apart from the afore- 
mentioned difference between the time-dependent and the sta- 
tionary models, a dependence of Afi on Rm cannot be inferred. 

In addition, the growth rate of the leading eigenmode of D, 
A - %icr), may be taken as a measure for the accuracy of the 
mean-field description. Ideally, it is 0, whereas all overtones are 
highly diffusive ( Schrinner et al. 2010a, 2011b). For numerical 
simulations, however, it is impossible to hit the critical point ex- 
actly. The growth rates of the fundamental mode for all mod- 



50 100 150 
Rm 

Fig. 2. Afi versus the magnetic Reynolds number for the (quasi)- 
stationary benchmark dynamo (triangle) and the 6 chaotically 
time-dependent dynamos considered. 



els are given in Table [T] in units of rj/L^. Note that the turbu- 
lent diffusivity largely exceeds the molecular one. For model 2, 
for instance, we find values of jS-components up to 63//. Thus, 
1/A ^ ip L^/t] for model 2 is in fact much larger than one ef- 
fective decay time and the fundamental mode is already near its 
critical state. With decreasing Afi, the growth rates become even 
smaller and the values closest to zero have been obtained for 
models 5, 6 and 7. 

For model 7, the best parameterization was obtained and 
Afi = 1.1 was achieved. This still is clearly larger than for the 
simulation of r otating magneto convection for which Afi = 0.28 
was found by ISchrinner et al.l (|2007). Figure |3] displays £ in 
comparison to Differences are visible in all components, 
but also principal similarities. 

Figure |4|compares the azimuthally and time-averaged mag- 
netic field resulting from direct numerical simulations with the 
fundamental eigenmode of D for model 7. Apart from small dif- 
ferences, direct numerical simulations and mean-field calcula- 
tions agree very well for this example. 
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Fig. 3. Contour plots of Sr, &e, and 6^ (top line, from left to 



right) and £*'\ £g *, and fi^'' (bottom line, from left to right) for 

model 7 in units of (v/L) (^JjUt/Q)^^^. In these units, contour lines 
correspond to +3, +9, +15. 
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6. Discussion 

We do not derive mean-field coefficients for time-averaged mean 
fields in a formal sense. Instead, the comparisons in Fig. |3]and 
Fig.|4]are motivated by assumptions i) and ii) in Section[3] They 
are reasonably well fulfilled for the simulations considered here. 
For model 7, v time-averaged over 30 magnetic diffusion times 
drops to 0.1% of V. A similar decrease has been found for b. 
Moreover, the axisymmetric flow and the axisymmetric mag- 
netic field show comparatively little time variation. The time- 
dependent components of V and B are of the order of 8% (37%) 
and 17% (23%) of the total (axisymmetric) flow and the to- 
tal (axisymmetric) magnetic field, respectively. In other words, 
time-averaging results in axisymmetric fields and azimuthal av- 
eraging reduces somewhat the time variability. Combining both 
leads to an extended averaging and thus to an improved scale 
separation. This finding is also supported by Fig. |5] and Fig. 
|6] Figure |5] shows a noticeable decrease of A£ for model 5, if 
the time-averaging period t is increased. A£ drops rapidly un- 
til T ^ 2L^ jri; an extension of the time-averaging interval over 
more than 2 magnetic diffusion times improves the scale sep- 
aration only slightly. Moreover, figure |6] compares the resid- 
ual with the mean magnetic field for model 1 (upper panel) 
and model 7 (lower panel). In contrast to model 1, the time- 
dependent residual and the time-averaged mean magnetic field 
of model 7 vary on fairly different length scales. Thus, a bet- 



Fig. 4. Comparison between direct numerical simulations and 
mean-field calculations for model 7. Bottom line: azimuthally 
and time-averaged magnetic field obtained by solving dglfTTTi. 
Top line: fundamental eigenmode resulting from (fT5t . Each 
component has been normalised with respect to its abso- 
lute maximum. Therefore, the greyscale coding varies from 
-1, white, to +1, black, and the contour lines correspond to 
+0.1, +0.3, +0.5 ±0.7, +0.9. 



ter scale separation for the time-dependent models seems to be 
plausible. Consequently, A£ drops drastically from model 1 to 
model 2. However, there are still noticeable differences between 
£ and £*^* in Fig.[3]for model 7. From this comparison alone, it 
would be difficult to decide, whether the parameterization of £ 
is already satisfactory. Further support for its reliability comes 
from the good agreement between direct numerical simulations 
and mean-field modeling in Fig. |4] Contour plots of all three 
components of the leading eigenmode of D agree well with cor- 
responding results from direct numerical simulations. Moreover, 
the growth rates of the leading eigenmodes are close to zero, as 
it has to be expected for stochastically stationary dynamos in the 
kinematically stable regime. 

Dynamo action in stars and planetary cores most probably 
takes place on a wide range of length and time scales. Due to 
computational limitations, numerical simulations do not cover 
the whole range of scales possibly involved. Global dynamo 
simulations focus on large scales only, and small-scale dynamo 
action as reported for box simulations (e.g. IVogler & Schiissleri 
2007) is typically not present. Therefore, the component of the 
electromotive force independent of the mean magnetic field, i.e. 
£*"^ in dUi, is zero for most global dynamo simulations. The 
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Fig. 5. A£ for model 5 varying with increasing time-averaging 
period t. 




Fig. 6. Upper panel: non-axisymmetric radial magnetic field at 
arbitrary longitude (left) and azimuthally-averaged radial mag- 
netic field (right) for model 1 . Lower panel: non-axisymmetric 
and time-dependent radial magnetic field at arbitrary longitude 
and time (left) and azimuthally and time-averaged radial mag- 
netic field (right) for model 7. The style of the contour plots is 
the same as in Fig.|4] 

possibility of a contribution to fi from small-scale dynamo ac- 
tion has ijeen intensively discussed in the litera ture (e.g. ,Radle3 
ll976HRadlei<l2000HRadler & Rheinhardl2007h . Recently, it has 
been used to argue against the applicability of mean-field the- 
ory (Cattaneo & Hughes 2009). However, conceptual difficulties 
which might result from simultaneous small- and large-scale dy- 
namo action are at present not relevant for the approach followed 
here. 

The problem of scale separation has been addressed in this 
study with regard to spatial scales, only. This seems to be ap- 
propriate because we intend to describe stochastically stationary 
mean-fields with zero growth rate resulting from self-consi stent 



dynamo simulations. More generally, time derivatives of the 
mean field in the expansion of £ have to be taken into ac- 
count. Similarly to non-local effects (Brandenburg et al. 2008]), 
also memory effects of the flow may influe nce the pa rameteriza- 
tion of the electromotiv e force dHubbard & Brandenburgw2009t 
iHughes & Proctodl2010l) . 

Mean-field theory as presented here is intrinsically a kine- 
matic approach. In general, an eigenvalue calculation based 
on ( fTsT i will lead to growing modes, even if the mean- 
field coefficients are derived f rom a saturated ve l ocity field 
(iTilgner & Brandenburg l2008l ICattaneo & TobiasI l2009b . In 
such a situation, mean-field models might not match direct 
numerical simulations, unless a more self-consistent exten- 
sion of ffi e ffieoretical framework and the test-field method 
is app lied (ICourvoisier et al ] 120 lot iRheinhardt & Brandenburg! 
l2010h . However, the class of chaotically time-dependent dy- 
namos we considered is kinematic ally stable. iSchrinner et al] 
(I2010ah identified a regime of fast rotating dynamos which 
are dominated by only one dipolar mode at marginal stability, 
whereas all overtones are highly diffiisive. A saturated velocity 
field from this class of dynamos does not lead to exponential 
growth of the magnetic field in a corresponding kinematic cal- 
culation. Consequently, the mean-field approach based on (fTSl l 
is applicable for this dynamo regime, as confirmed once more 
in this study. In particular for models beyond this regime, the 
reliability of the mean-field approach presented here is not guar- 
anteed and has to be made plausible by a comparis on with direct 
numerical simulations (e.g. lSchrinner et ani2011al) . 

7. Summary 

The mean-field description of a (quasi-)stationary dynamo 
suffered from a poor scale-separation and therefore from 
an insufficient param eterization of the electromotive force 
(ISchrinner et al.ll2b07l) . For the chaotically time-dependent dy- 
namos considered here, both improves a lot, if a combined az- 
imuthal and time average is applied. The more accurate param- 
eterization of £ leads to a good agreement between mean-field 
models and direct numerical simulations: Field topologies and 
growth rates resulting from both approaches are very similar. 

In conclusion, the test-field method to determine mean-field 
coefficients from direct numerical simulations proves to be re- 
liable for chaotically time-dependent dynamos, too. Mean-field 
theory may serve as a powerful tool to analyse dynamo processes 
in global models resulting from direct numerical simulations. 
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